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Abstract 

We study the quantum decay of DO-branes in two-dimensional OB string 
theory. The quantum nature of the branes provides a natural cut-off for the 
closed string emission rate. We find exact quantum mechanical wavefunctions 
for the decaying branes and show how one can include the effects of the Fermi 
sea for any string coupling (Fermi energy). 
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1 Introduction 



Originally dynamical triangulations of string world-sheets were introduced as 
a reparametrization invariant lattice regularization. It successfully defined 
non-critical strings of fixed world-sheet topology in a non-perturbative way, 
the scaling to the continuum controlled by the lattice spacing PUEHS]- The so- 
called matrix model description was a very convenient way of implementing 
the combinatorial task of summing over all abstract triangulations of the 
world-sheet, which automatically, via the large N expansion, arranged the 
world-sheets according to topology. No physical interpretation of the matrix 
itself was given until recently when it was suggested in the case of the c = 1 
matrix model that the matrix could be given the interpretation as the open 
string tachyon field between N unstable D-branes (N being the size of the 
matrix), the unstable D-branes themselves being identified as the eigenvalues 
of the matrices |U E] . 

This intriguing picture has passed a number of non-trivial tests (also for 
models with c < 1, [HI 13 El) and it offers the possibility for the first time to 
study quantum D-branes in strongly coupled string theories. 

In this paper we will discuss the decay of such quantum D-branes. 

The outline of this paper is as follows. In section 2 we describe the decay 
of the DO-brane in the classical approximation and reproduce directly from 
the classical motion the closed string tachyon 1-point function. In section 3 
we present the exact quantum mechanical treatment of the same process, 
although without taking into account the effect of the Fermi sea, and show 
how a natural cut-off arises for the closed string emission. In section 4 we 
show how to exactly incorporate the effects of the string interactions (Fermi 
sea) in the preceeding quantum mechanical picture. We close the paper with 
a discussion. 

2 Classical decay of the DO-brane 

In the double scaling limit the ground state of 2d bosonic string theory is con- 
structed by filling one side of the upside-down harmonic oscillator potential 
—X 2 /2a' with free fermions up to a Fermi level —fi, where zero energy is the 
top of the potential. The string coupling g s ~ and the partition func- 
tion, closed string tachyon operators, macroscopic loop operators etc have a 
unique perturbative expansion in g s which can be obtained from the exactly 
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solvable quantum mechanics of the upside-down harmonic oscillator. As first 
pointed out in j^j the situation is unclear if we move to strong string cou- 
plings, i.e. small /i near the top of the potential. Clearly tunneling (so-called 
"non-perturbative" effects) between the two sides is no longer exponentially 
suppressed, but there is no "first principle" in the bosonic case telling us how 
to relate the two sides. 

Fortunately in the case of two-dimensional OB and OA superstrings, as 
realized in (TUJ HI] , this ambiguity is lifted and both sides are filled up to the 
same level. 

In the following we will perform calculations within the OB model 3 defined 

by 

JdTe-f dt {^ 2+v ^} (1) 

where the nondynamical gauge field in the covariant derivative just restricts 
the path integral to singlet states and so the standard free fermions give a 
complete description of all the degrees of freedom. 

Further, the quantum states of a DO-brane are precisely the quantum 
states of the Hamiltonian of the inverted harmonic potential except that the 
spectrum starts at — fi and the model provides us with a complete description 
of the dynamics of a single DO-brane |13j . 

The potential for the OB matrix model differs from the pure bosonic 
potential by a factor of two: 

V = — —A 2 (2) 

As realized in |3] eigenvalues may be identified with (unstable) DO-branes. 
Consequently time dependent solutions are of interest since they may be iden- 
tified with decay processes. A solution of the classical equations of motion 
is 

t 



A(t) = y 4/io;' sin 7rA cosh — ^==, (3) 

where A should not be confused with the eigenvalue A. Its energy relative to 
the Fermi surface Ep = —\l is /i cos 2 7rA. This solution is the matrix model 
analog of Sen's 'rolling tachyon' ^3] and it will emit closed strings tachyons, 
the rate dictated by its contribution to the on-shell closed string 1-point 

3 However, all results derived here away from the strong coupling region (g s large) are 
qualitatively correct also for the the 2d bosonic string (the c = 1 matrix model) . 
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function. This in turn can be computed since the closed string tachyon (in 
the NS-NS sector) can be defined ^U] in the OB matrix model as 4 

T NS -ns(E) ~ lira (leg factor) • / dt e iEt (tr e~ m ^ ) (4) 

where the so-called leg-factor is 

leg factor = , (5) 

T(-iEJa'/2) 



and the expectation value is with respect to the (double scaling limit) matrix 
integral: 

(tre-^ 2 W)= fd\p ds (X,t)e- ilx2 , (6) 



Pdsi^yt) being the appropriate double scaling limit of the eigenvalues. 

Let us compute the contribution of the classical motion A(t) of the eigen- 
value (j3J) to the 1-point function of the NS-NS closed string tachyon. The 
contribution of such a classical eigenvalue to p<fe(A,t) is simply 5(\ — \(t)). 
Thus we have to evaluate the integral 

dt e iEt e ~ mm ' sin2 wX cosh2 vh (7) 

in the limit I —> 0. The non-trivial small I behavior comes from the the large 
t region of the integral and the following change of integration variable 

~2~ 



u = leV^ (8) 
leads to an integral 

l ^Lf iE ^i J™ ^g-Wusm^y^y ~ (9) 

Performing the integral, setting a' = 2 and inserting the leg-factor (j3J) we 
obtain the correct answer: 

-»Elogsin 2 7rA -iE 

V{-iEf ' { ' 

Analogous calculations here and below, can of course be easily done also 
for the closed string tachyons in the R-R sector. 

4 In we have made a rotation I —> il, in accordance with the discussion in |15j . 
sec. 11. 
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3 'Quantum' decay - D-brane wave packets 



As already mentioned the quantum theory of a DO-brane in the c = 1 theory 
is completely described as a free fermion with the Hamiltonian 

Id 2 1 

H = —^---K 2 \ 2 (11) 

where k = 1/ s/la' . Therefore the classical picture of DO-brane decay is only 
approximate, limited by the uncertainty principle. In the quantum theory we 
are forced to consider wave packets instead of localized eigenvalues following 
a classical trajectory with both a definite position and momentum. The 
only non-trivial aspect is the fermionic nature of the DO-branes, which gives 
an indirect interaction with the Fermi-sea background and thus reflects the 
nonzero value of the string coupling. 

Let ip(\,t) be a normalized solution to the Schrodinger equation for H 
and let us assume it has no overlap with the eigenfunctions of H with energy 
E < Ep- The contribution of this quantum state to the density Pd s (\ t) will 
be \ip(\,t)\ 2 and thus the corresponding contribution to when inserted 
in (J3J), leads to the following integral: 

T q Ns-Ns( E ) = j im (leg factor) • J dte iEt J d\e~ ilx2 \^(X,t)\ 2 . (12) 
For an arbitrary solution to the Schrodinger equation we can always write 

^(\,t) = J2c E M\t). (13) 

E 

If this expansion of if) contains E < Ep they have to be cut away and the 
wave function re- normalized (see section 4 below). 

Let us concentrate on the simplest wave packets, Gaussian wave packets, 
and ignore at first the problem of overlap with the Fermi sea. We will address 
it in the next section. First a few general observations: (1) for a quadratic 
potential the expectation values (A(t)) will always follow a classical orbit 
as follows from Ehrenfest's theorem. (2) A wave packet which is Gaussian 
at some time t will remain Gaussian. This follows because the propagator 
G(\, A'; t) corresponding to H is Gaussian in A, A'. (3) Since the peak will 
then coincide with the expectation value of A we know that for an initial 
wave packet of the form 

^,(A, 0) = {^j * e ~-(^of+ iP0 x ( 14 ) 
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we have 



MA,f)| 2 = e -Mt)(x-x(t)? (is) 



7T 



where X(t) is just the classical orbit corresponding to initial values of A ,Po 
and a(t) can be calculated to be 5 

4 ^ 

a(t) = — — , A(t) = cosh 2 nt + — sinh 2 nt (18) 

ZA ( 6 J 

Consequently the temporal evolution of wave packets is essentially dictated 
by the classical energy 6 i7 c ;(A ,Po) — f G°o — k2 ^o) since the peak just follows 
the classical orbit and the wave packet never splits in two, as a generic wave 
packet would do in the inverted harmonic potential. However, due to the 
very rapid spread of the wave one still has non-zero transmission though the 
potential barrier even if H c i < 0. For a given initial wave packet we see that 
a true classical picture of a localized wave function is only valid for times 
times Kt < log a or t < y/c/ log a. 

Let us for simplicity consider a wave packet where po = and Ao > 0. 
This is the wave packet version of Sen's rolling tachyon and Ao will now be 
related to \i as follows: 

Aq = 4/ic/ sin 2 7rA. (19) 

We can now calculate the closed string tachyon one-point function (|12|) . The 
integral over the eigenvalues gives in this case 

1 / , A 2 , cosh 2 , . 

exp ( - % l \ » ) (20) 
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It might seem as if one could neglect the contribution I A /2a from the spread- 
ing of the wave packet since we take the / — > limit. However since A grows 



5 It can easily be read off from the propagator 



G(A, A'; t) = ^- V2 exp ( , ^ [(A 2 + A' 2 ) cosh Kt - 2AA'l S> , (16) 



sinh nt 



ip{X,t) = J dX' G(A,A';t)V(A',0). (17) 
3 The actual quantum energy of the wave packet is (ip\H\ip) = H c i(X ,p ) + f — fr- 
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exponentially in time this contribution is of the same order as the classi- 
cal piece. Therefore the classical result (|1U|) will get modified. Making the 
change of variable (JBJ) we have to perform now the integral 



X 2 

i Ey /a'/2 f°° du lE J^T2 . e-^* 



a'/2l-^ a ' 2 / — u l ^ a ' 2 ■ - (21) 

Jo U y'l + iCu 



where 

1/1 2a\ 



C = l(L + ^) = -J^m l — + 4aJrt/2\ . (22) 



v 4ayaV2 



This integral can be performed by a further change of variable v 
which leads to an integral of the form (with a' = 2) 



iCu 
1+iCu 



UC)- iE f 1 dvv iE 'Hl - v )- iE -h-^- v (23) 
Jo 



The final result, after including the leg-factors, is 
quantum - 



F(iE) 




F(-iE) 1 




T(iE) 


If, 


T(-iE) 1 





- iE '1 .„\ 1 „ (.„ 1 



2 / V 2 

classical = j^r^j ■ I y J (25) 

where for comparison we have also written the classical result and where 
the relation between Ao and /i is as in ()19j) . Note that this appearance of 
// (equivalently g s ) here is purely 'kinematical' and enters only through the 
initial condition Ao- The appearance of \x due to string interactions (influence 
of Fermi sea) will be treated in the next section. 
Let us first discuss the case when A 3> 0. 

First note that from the asymptotics \F\ (a, 6; — x) ~ x~ a Y{b) /T(b — a) 
we see that the quantum result for the tachyon 1-point function leads to the 
classical one in the limit \i — > oo (g s — > 0) while keeping the energies E fixed. 

However for any large but finite \x the behavior of both formulas is quite 
different. The classical result is a pure phase factor which for large E behaves 
as: 

Tns-ns(E) ~ • -)~ iE ■ e - 2iE e 2iElogE + 0(1/ E) (26) 

Thus the number of emitted particles, N = J ^-\T NS -ns(E)\ 2 , diverges loga- 
rithmically and the expectation value of the emitted energy diverges linearly 
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[5] as in the similar calculation in 26 dimensions [16]. It was suggested that 
a cut-off of order l/g s should be put into these calculations but in our case 
of quantum DO-branes we see that such a cut-off arises naturally. Indeed, 
denoting M = Aq/(4C) we have the large E asymptotics 

(iE, ~; -M^j ~ e~% cosh ^VM^-iE^j (27) 
V(l/2 -iE) ~ V2^e-^ E e iE - iElogE (28) 



We see that the answer is cut off when tcE/2 ^> y / 2ME, or, dropping con- 
stants, when E 3> M. Thus the energy emitted is finite, the regularization 
provided by the quantum mechanical nature of the DO-brane as indeed conjec- 
tured in j^j. If the classical orbit has a turning point away from zero (Ao 3> 0) 
then the cut-off is (with a' reinserted and \i\fa' ~ l/<7«) 

yfc'Ecaof, ~ - = l - =- ~ -, (29) 

caya ' + (1/caya') 9s 

where the constant c is of order 1 and the last ~ is valid when the location 
of the wave packet at t = is of the order of \fa'. Anyway, l/g s will always 
serve as a upper cut-off as long as the A 3> 0. 

Let us now turn to the opposite case. When the turning point of the 
classical orbit is close to the maximum of the potential M goes to zero and 
the cut-off is now set just by ([28)1 . Then in the r.h.s. of (J2l?|) l/g s should 
be replaced by 1, showing that quantum effects wash out the signature of 
string perturbation expansion. In particular this is the case in the "pure 
quantum case" where Ao = and the classical eigenvalue is located on the 
top of the potential without rolling down. This situation is not sustainable 
for the quantum brane and the probability of finding the eigenvalue within 
a distance d from the origin decreases like e~ Kt d\fa for the Gaussian wave 
packet. The cut-off of the emitted energy is then just E cuto ff = \j\ja! . 



4 Large g s and inclusion of the Fermi sea 

Above we obtained the exact motion of the wave packet in the inverted 
potential. However we treated the motion independently from the other 
eigenvalues. This is a good approximation if the overlap with the Fermi sea 
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is small. If we consider wave packets imitating to some degree the rolling 
tachyon of Sen we expect the above calculation to be reliable if iiyfa' ~ l/g s 
is large. However, when the string coupling is large the Fermi level is close 
to the top and we are bound to get a significant overlap between the wave 
packet and the Fermi sea and it has to be taken into account. 

The exact iV-body wave function of the system will just be a Slater de- 
terminant with the ip(\,t) wave function derived in the previous section 
as one of its components. Then any observable such as the closed string 
tachyon 1-point function has to be computed using the Slater determinant 
wave function. Such a calculation would be in general quite formidable due 
to the needed anti-symmetrization and lots of possible subtle cancellations. 
This will introduce another source for the dependence on the string coupling 
g s into the picture (since this is encoded in the Fermi level). 

In order to bypass these complications we will construct from if)(\,t) 
an equivalent /z-dependent wave function ip e ff(\,t) which will render anti- 
symmetrization trivial and allow to use single particle intuitions. 

To this end let us consider properly normalized Slater determinants. If all 
the component wave functions are orthogonal to each other then the Slater 
determinant will be properly normalized. This is certainly the case for the 
levels in the Fermi sea, but the wave packet ip(X,t) will have components 
also below the Fermi level. Of course these components will not contribute 
to the iV-body wave function. So the only effect of the Fermi sea will be to 
truncate the original expansion of the wave packet 

V(A,t) =J2c E e- tEt M\) (30) 

E 

to 

^ eff (X,t) = J\f ■ ?p proj (\,t) —M £ c E e- iEt 4> E {\) (31) 

E>-n 

where the normalization constant is 

A/-= [ £ l^l 2 ) 2 (32) 

\E>-n J 

Note that the normalization constant M will re-normalize all expectation 
values (average energies, 1-point functions) from the single particle case. 
The energy profiles of the 1-point functions will of course also change due to 
the absence in ip pro j(\, t) of some of the original components of ^(A, t). Since 
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Norm E 




Figure 1: The normalization constant M and the energy of the effective 
wave-function if) e ff(\,t) (corresponding to the Gaussian wave packet (f55|)) 
as a function of //. 



ip e ff(\, t) is orthogonal to all states below the Fermi level anti-symmetrization 
is trivial and effectively drops out. 

The projector is Pe f (E) = 6(E — Ep), the Fourier transform of which is 

1 j p iE F {t-s) 

P^s) = -S { t- s) + ^ r - ) . (33) 

Thus we can write (with Ep = —ji) 

I ? r°° d s r i 

^proj(\t) = -ip(\,t) + — - \e-^{X,t + s)-e^(X,t-s) . (34) 



2 rv ' 2vrJo s 

This integral representation is of course completely general and might be 
convenient whenever one actually knows the wave function ip(\,t). 

The simplest case, and the one in the family of states we have considered 
here which is least semi-classical, is 



I A 2 lie 4 1+itanh 2 

MA,0) = — -re-T, Vo(A,t) = —__==—_=- (35) 
(2?r) 2 (2tt) i ,/cosh | ,/l + 1 tanh f 



where the last expression follows from ()17|) . With this choice of wave func- 
tion both the classical energy H d = \{j>q — k 2 Aq) and the quantum energy 
(^ |##o) = H d (\ ,p ) + f - are zero. 

In order to extract ip e ff(\, t) from t/>o(A, £) one can use the general formula 
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mu=l . 5 mu=l 




Figure 2: The wave-functions i[) e ff(\,t = 0) (dotted points) corresponding 
to the Gaussian wave-packet i[>o(\,t = 0) in eq. (|55j) (shown as a solid line) 
calculated for \i = 1.5, 1.0, 0.5 and 0. 



3 Inn / . \ 3 



,2 



to get after some manipulations 7 : 
- e'^^e^ 

v^tt! J -i t + zf + log^ 

(36) 

Alternatively in this particular case one can explicitly find the decomposition 
in energy eigenstates 8 : 



r(i/4- 


-iE) 


r(3/4- 


-iE) 



^o(A,t) = / dEc E M^~ iM , c E = 2-^-—— (37) 

(47rcosh(27rE))4 

where the coefficients fall off exponentially for large \E\, ce ~ |i?| _1//4 e _7r ' £ '' //2 . 
For large /x the overlap with the Fermi sea will indeed be exponentially small 

7 Write l)34|) as f°° dfJ,-j-ip pro j , interchange the /i and s integrations and change variables 
from s to v. 

8 The parity even parabolic cylinder functions iPe(A) can be found in |12l I15j . we use 
here the normalization from ^H] and E = — | where a is the parameter used in |15| . 
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in fi and thus "non-perturbative" in nature. However, the "non-perturbative" 
nature is of the same origin as most other "non-perturbative" corrections 
discussed in the literature since the exponential nature of the correction 
comes from the Fermi-sea wave functions being in a classically forbidden 
region (around A = 0). The same can be said about the normalization 
constant M from (J32j) . Note also that the energy of this "quantum" brane 
will change from zero to a positive value due to the interaction with the Fermi 
sea. 

In Fig. ^ we have shown the behavior of the normalization constant M 
and the energy E e ff defined by 

Eeff = ^eff\H\i> eff ) . (38) 

Only for ^/a'n < 1, i.e. for g s > 1 is there an effect which is not exponentially 
suppressed in \i. This is corroborated by looking at the wave function ip e ff(\) 
itself. However, for small fi we see quite a large change in the wave function 
due to the interaction with the Fermi shown in Fig. 

5 Discussion 

This work was motivated by the calculations of closed string emission in the 
time-dependent rolling tachyon open string background of Sen. These calcu- 
lations, in the context of 26-dimensional bosonic string theory were performed 
in and in the context of 2d critical string theory and/or matrix models 
in El El HOI • The observation was that the emitted energy was generically 
infinite. In the concrete calculations the open string background was viewed 
as purely classical 9 , which in the matrix model formulation translates into 
the statement that the eigenvalue corresponding to the DO-brane follows a 
classical trajectory in the inverted harmonic potential. Indeed, we saw that 
it was very simple to perform the closed string emission rate calculation using 
this classical trajectory. Here again the emitted energy is infinite, in accor- 
dance with the fact that the DO-brane is treated as a classical background 
object. 

Since the quantum states of a DO-brane are in a one to one correspon- 
dence with the quantum states of the cut-off Hamiltonian in the inverted 

9 In |S] there is a discussion of using bosonization of chiral fermions to get a more 
complete treatment. 
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harmonic potential, the matrix model offers a simple way to move away from 
the classical situation. In this paper we derived an exact quantum descrip- 
tion of the time-dependent decay process of the unstable DO-brane in type OB 
two-dimensional superstring theory. By doing so the emitted energy indeed 
becomes finite, and for small g s it is cut off at \fa'E = l/g s . For large g s the 
cut-off is of order 1. 

As often before, we are in a situation where the matrix model offers 
us a simple way of addressing questions which are not easily addressed in 
the continuum or higher dimensional theory. In particular here we could 
exactly treat a quantum open string background since, as discussed above, 
the quantum mechanics of fermions in the upside-down harmonic potential 
with energy levels below the Fermi level at Ep = —fj, filled should be viewed 
as a candidate for the continuum quantum open string field theory in 2d, 
describing the dynamics of DO-branes. 

In order to apply the lesson from matrix models to higher dimensional 
critical string theories one needs to understand how to describe the quantum 
D-brane (quantum open string background) from a string point of view. One 
possibility is to use Witten's open string field theory at the quantum level. 
However it would be most attractive to have a direct worldsheet description. 
Presently we can understand that the classical D-brane should be associated 
with a boundary conformal field theory, but the quantum D-brane generalizes 
this situation, and there should be a suitable string-theoretical description of 
the quantum D-brane. The ease with which the situation is handled in the 
matrix model context and the fact that very sensible results emerge is a strong 
hint that there should exist a simple higher dimensional string description 
too. 
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